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Abstract 

We report the measurement of sin 2 9yy in v — N deep inelastic 
scattering from the NuTeV experiment. Using separate neutrino and 
anti-neutrino beams, NuTeV is able to determine sin 2 9y/ with low 
systematic errors by measuring the Paschos-Wolfenstein variable R~, 
a ratio of differences of neutrino and anti-neutrino neutral-current and 
charged-current cross-sections. NuTeV measures sin 2 8w( on ~ shcl1 ^ = 
0.2253 ± 0.0019(stat) ± 0.0010(syst), which implies M w = 80.26 ± 
0.11 GeV. 



1 Introduction 



In the past, neutrino scattering experiments have played a key role in 
establishing the validity of the electroweak Standard Model. Today, even 
with the large samples of on-shell W and Z bosons at e + e~ and pp colliders, 
precision measurements of neutrino-nucleon scattering still play an impor- 
tant role. Not only are these measurements competitive in precision with 
direct probes of weak boson parameters, but they also test the validity of the 
electroweak theory in different processes and over many orders of magnitude 
in q 2 . In this respect, if neutrino scattering observed deviations from expec- 
tations based on direct measurements from W and Z bosons, this would be 
an exciting hint of new physics entering in tree-level processes or in radiative 
corrections. In particular, neutrino scattering would be sensitive to a differ- 
ent menu of non-Standard Model effects ranging from leptoquark exchange 
to neutrino oscillations [||, 

Experimental quantities sensitive to electroweak physics that are most 
precisely measured in neutrino scattering are the ratios of charged-current 
(W exchange) to neutral-current (Z exchange) scattering cross-sections from 
quarks in heavy nuclei. The ratio of these cross-sections for either neutrino 
or anti-neutrino scattering from isoscalar targets of u and d quarks can be 
written as|| 

<t(V„ JV -. ti-MX) 



where 



a^N -> fi+X) 1 



and g\ R = u\ R + d\ R , the isoscalar sums of the squared left or right- 
handed quark couplings to the Z. At tree level in the Standard Model, 
q L = I^eak~ Qemshi 2 6\y and g# = — Qemsui 2 9 w ; therefore, R v is particularly 
sensitive to sin 2 #jy- 

In a real target, there are corrections to Eqn. [I] resulting from the presence 
of heavy quarks in the sea, the production of heavy quarks in the target, non 
leading-order quark-parton model terms in the cross-section, electromagnetic 
radiative corrections and any isovector component of the light quarks in the 



target. In particular, in the case where a charm-quark is produced from scat- 
tering off of low-x sea quarks, the uncertainties resulting from the effective 
mass suppression of the heavy final-state charm quark are large. The un- 
certainty in this suppression ultimately limited the precision of previous vN 
scattering experiments which measured electroweak parameters^, [5|, ||. 

To eliminate the effect of uncertainties resulting from scattering from sea 
quarks, one can instead form a quantity suggested by Paschos and WolfensteinfTj, 



Since o~ vq = o~ vq and a vq = o~ uq , the effect of scattering from sea quarks, which 
is symmetric under q <-> q, cancels in the difference of neutrino and anti- 
neutrino cross-sections. While allowing substantially reduced uncertainties, 
R~ is a more difficult quantity to measure than R" , primarily because neutral 
current neutrino and anti-neutrino scattering have identical observed final 
states and can only be separated by a priori knowledge of the initial state 
neutrino. 

2 The NuTeV Experiment and Neutrino Beam 

The NuTeV detector consists of an 18 m long, 690 ton target calorimeter 
with a mean density of 4.2 g/cm 3 , followed by an iron toroid spectrome- 
ter. The target calorimeter consists of 168 iron plates, 3m x 3m x 5.1cm 
each. The active elements are liquid scintillation counters spaced every two 
plates and drift chambers spaced every four plates. There are a total of 84 
scintillation counters and 42 drift chambers in the target. The toroid spec- 
trometer is not directly used in this analysis. NuTeV used a continuous test 
beam of hadrons, muons and electrons to calibrate the calorimeter and toroid 
response. 

In this detector v^jv^ charged-current events are identified by the pres- 
ence of an energetic muon in the final state which travels a long distance in the 
target calorimeter. Quantitatively, a length is measured for each event based 
on the number of neighboring scintillation counters above a low threshold. 
Charged-current candidates are those events with a length of greater than 20 
counters (2.1 m of steel-equivalent), and all other events are neutral-current 
candidates. 

NuTeV's target calorimeter sits in the Sign-Selected Quadrupole Train 
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(SSQT) neutrino beam at the FNAL TeVatron. The observed neutrinos 
result from decays of pions and kaons produced from the interactions of 
800 GeV protons in a production target. Immediately downstream of the 
target, a dipole magnet with / Bdl = 5.2 T-m bends pions and kaons of 
one charge in the direction of the NuTeV detector, while oppositely charged 
and neutral mesons are stopped in dumps. Focusing magnets then direct the 
sign-selected mesons into a 0.5 km decay region which ends 0.9 km upstream 
of the NuTeV detector. The resulting beam is either almost purely neutrino 
or anti-neutrino, depending of the selected sign of mesons. Anti-particle 
backgrounds are observed at a level of less than 1-2 parts in 10 3 . The beam 
is almost entirely muon neutrinos, with electron neutrinos creating 1.3% and 
1.1% of the observed interactions from the neutrino and anti-neutrino beams, 
respectively. 

Because charged-current electron neutrino interactions usually lack an en- 
ergetic muon in the final state, they are almost always identified as neutral- 
current interactions in the NuTeV detector. Therefore, the electron neutrino 
content of the beam must be very precisely known. Most (93% in the neu- 
trino beam and 70% in the anti-neutrino beam) observed v e /v e s result from 
decays. The remainder are products of prompt decays of charmed par- 
ticles or neutral kaons, or decays of secondary muons. Prediction of the 
former component comes from a beam Monte Carlo, tuned to reproduce the 
observed v^jv^ flux (Figure [I]). Because of the precise alignment of the mag- 
netic optics in the SSQT, this procedure results in a fractional uncertainty on 
the prediction of v e jv e from of ~ 1.5%, dominated by the K% branch- 
ing ratio uncertainty. Small detector calibration uncertainties, 0.5% on the 
calorimeter and muon toroid energy scale, affect the measured v^jv^ flux 
and also contribute substantial uncertainties to both the muon and electron 
neutrino fluxes. Sources of v e jv e other than K ± decay have larger uncertain- 
ties, at the 10-20% level, because of the lack of a direct constraint from the 
data. 

3 Extraction of sin 2 6w 

Events selected for this analysis are required to deposit at least 20 GeV 
in the target calorimeter to ensure efficient triggering and vertex identifica- 
tion. The location of the neutrino interaction must be within the central 
2/3 rds of the calorimeter's transverse dimensions, at least 0.4 m of steel- 
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1: The and z/^ energy spectra from the data and the tuned beam 
Carlo. 



1.1 Data/MC 
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Figure 2: Length distributions in the data from the neutrino and anti- 
neutrino beams. Neutral-current /charged-current separation is made at a 
length of 20 counters, approximately 2.1 m of steel. 



equivalent from the upstream end of the calorimeter, and at least 2.4 m from 
the downstream end. The first requirement reduces the misidentification of 
^mAv even ts with muons exiting the side of the calorimeter; the second re- 
duces non-neutrino backgrounds, and the third ensures sufficient calorimeter 
downstream of the interaction to measure the event length. Small back- 
grounds from cosmic-ray and muon induced events are subtracted from the 
sample. After all cuts, 1.3 million and 0.30 million events are observed in the 
neutrino and anti-neutrino beam, respectively. The ratios of neutral-current 
candidates (short events) to charged-current candidates (long events), R meas , 
are 0.4198 ± 0.0008 in the neutrino beam and 0.4215 ± 0.0017 in the anti- 
neutrino beam. 

R meas is related to the ratios of cross-sections and sin 2 Q w using a de- 
tailed detector and cross-section Monte Carlo simulation with the tuned flux 
(Figure Q) as input. This Monte Carlo must predict the substantial cross- 
talk between the samples. In the neutral-current sample, the backgrounds in 
the neutrino and anti-neutrino beam from v^jv^ charged-current events are 
19.3% and 7.4%, and the backgrounds from v e jv e charged- currents are 5.3% 
and 5.8%. The charged-current sample has only a 0.3% background from 
neutral-current events for each beam. 

The important details of the detector for this analysis are the calorimeter 
response to muons, the measurement of the neutrino interaction vertex, and 
the range of hadronic showers in the calorimeter. The efficiency, noise and 
active areas of the scintillation counters are all measured using neutrino data 
or muons from the testbeam. Longitudinal and transverse vertex resolutions 
and biases are studied using a GEANT-based detector Monte Carlo. The 
longitudinal bias, which directly impacts the length measurement, is mea- 
sured from the data using track-based vertices in events with two energetic 
final state muons. Hadronic shower length in the calorimeter is measured us- 
ing hadrons from the testbeam. To study possible effects from the difference 
in strange-quark content between neutrino-induced and 7r~-induced showers, 
hadronic showers from K~s are used as a cross-check. No significant differ- 
ences are observed. Measured detector parameters are varied within their 
uncertainties in the Monte Carlo to study systematic errors associated with 
this simulation. 

The cross-section model is of paramount importance to this analysis. 
Neutrino-quark deep- inelastic scattering processes are simulated using a leading- 
order cross-section model. Neutrino-electron scattering and quasi-elastic 



scattering are also included. Leading-order parton momentum distributions 
come from a modified Buras-Gaemers parameterization || of structure func- 
tion data from the CCFR experiment]^] which used the same target-calorimeter 
and cross-section model as NuTeV. The parton distributions are modified 
to produce u and d valence and sea quark asymmetries consistent with 



muon scattering ||10|| and Drell-YanfTl|| data. The shape and magnitude of 
the strange sea come from an analysis of events in CCFR with two oppo- 
sitely charged muons (e.g., vq — > fi~c, c — > /i + X)|L2|]. Mass suppression 
from heavy quark production is generated in a slow-rescaling model whose 
parameters are measured from the same dimuon data. The charm sea is 



taken from the CTEQ4L parton distribution functions ||13||. The magnitude 
of the charm sea is assigned a 100% uncertainty and the slow-rescaling mass 
for (y/v)c — > {yjv)c is varied from m c to 2m c . Our parameterization of 



-Riong = &l/vt is based on QCD predictions and datajnj and is varied by 



15% of itself in order to estimate uncertainties. Electroweak and pure QED 
radiative corrections to the scattering cross- sect ions are applied using com- 



puter code supplied by Bardin[l5], and uncertainties are estimated by varying 
parameters of these corrections. Possible higher-twist corrections are consid- 
ered with a 100% uncertainty using a VMD-based model which is constrained 
by lepto-production data[l6[ . 



The key test of the Monte Carlo is its ability to predict the length dis- 
tribution of events in the detector. Figure || shows good agreement between 
the data and Monte Carlo within the systematic uncertainties. 

To compute sin 2 9w, a linear combination of -R meas and -R meas was formed, 

-^meas = -^meas — ^^meas' (^) 

where a is calculated using the Monte Carlo such that -R~ eas is insensi- 
tive to small changes in the slow-rescaling parameters for charm produc- 
tion, a = 0.5136 for this measurement. This technique is similar to an 
explicit calculation of R~ , but here the background subtractions, the cross- 
section corrections to Eqn. ^, and the dependence on sin 2 9w are calculated 
by Monte Carlo. This approach explicitly minimizes uncertainties related to 
the suppression of charm production, largely eliminates uncertainties related 
to scattering from sea quarks, and reduces many of the detector uncertainties 
common to both the v and V samples. Uncertainties in this measurement of 
sin 2 9w are shown in Table [I]. 



SOURCE OF UNCERTAINTY 


5 sin 2 9w 


Statistics: Data 
Monte Carlo 


0.00188 

0.00028 


TOTAL STATISTICS 


0.00190 


Energy Measurement 
Event Length 


0.00045 
0.00051 
0.00036 


TOTAL EXP. SYST. 


0.00078 


Radiative Corrections 
Strange/ Charm Sea 
Charm Mass 
u/d, u/d 

Longitudinal Structure Function 
Higher Twist 


0.00051 
0.00036 
0.00009 
0.00027 
0.00004 
0.00011 


TOTAL PHYSICS MODEL 


0.00070 


TOTAL UNCERTAINTY 


0.0022 



Table 1: Uncertainties in sin 2 9w 
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Figure 3: Current direct My/ measurements compared with this result 



The preliminary result from the NuTeV data is f\ 



sin 2 9 w {on shell) = 0.2253 ± 0.0019(stat) ± O.OOlO(syst) 

/M top 2 - (175 GeV) 2 )\ 
-0.00142 x - top v J n 



(100 GeV) 



2 



+0 ,0048 x ,o ge ( . (5) 

The small residual dependence of our result on M top and Mmggs comes 
from the leading terms in the electroweak radiative corrections Since 
sin 2 e w ^ on - shcl ^ = 1 - M^/Mf, this result is equivalent to 

M w = 80.26 ±0.10(stat) ±0.05(syst) 

/ M top 2 -(175 GeV) 2 )\ 



+0.073 x 



(100 GeV) 



2 



A comparison of this result with direct measurements of is shown in 
Figure [3|. 

It is possible to extract the NuTeV result in a model-independent frame- 
work, where the result is expressed in terms of combinations of the left and 
right-handed quark couplings. The linearized constraint (expanded around 
one-loop couplings at an average log 10 ( \ ~^ e yi ) ~ 1 for NuTeV's central value 
of sin 2 6 W ) is 

0.4530 - sin 2 6 W = 0.2277 ± 0.0022 

= 0.8587*4 + 0.88284 - 1.1657w^ - 1.22884- ( 7 ) 

Note the similarity of this result to 1/2 — sin 2 9w = g\ ~ 9% the definition 
of the Paschos-Wolfenstein R~ in Eqn. [3[ 

*The weak radiative correction applied to extract sin 2 #yi/ on ~ shc11 from the measured 
quantities has changed since the presentation at Moriond due to an error in the implemen- 
tation of the Bardin code for radiative corrections. Two other small experimental correc- 
tions, for muon energy deposition and for charm semi-leptonic decays, were improved as 
well. The net shift in the result, 0.0054, is dominated by the fix in the implementation of 
the radiative corrections. 



(It is also possible to combine the NuTeV result with data from NuTeV's 
predecessor, the CCFR experiment. Adding the CCFR data@] in the -R^eas" 
based method described above, we obtain a slight improvement in precision, 
sin 2 6 W = 0.2255 ± 0.0018(stat) ± O.OOlO(syst).) 

4 Conclusions 

The NuTeV experiment has completed its data-taking and has extracted 
a preliminary result for sin 2 9y/ ( on_she11 ) , which is equivalent to My/ in the 
Standard Model. The precision of this result is approximately a factor of 
two improvement over previous measurements in z/N scattering because of 
the reduced systematics associated with measuring the Paschos-Wolfenstein 
ratio, R~ . This result is consistent with the average of direct My/ data. 

We would like to gratefully acknowledge the substantial contributions in 
the construction and operation of the NuTeV beamlines and the refurbishment 
of the NuTeV detector from the staff of the Fermilab Beams and Particle 
Physics Divisions. 

References 

[1] P. Langacker, et al. , Rev. Mod. Phys. 64, 87 (1991) . 

[2] K.S. McFarland, D. Naples, et al. , Phys. Rev. Lett. 75, 3993 (1995) . 

[3] C.H. Llewellyn Smith, Nucl. Phys. B228, 205 (1983) . 

[4] K.S. McFarland, et al. , Eur. Phys. Jour. CI, 509 (198) . 

[5] A, Blondel, et al. , Zeit. Phys. C45, 361 (1990) . 

[6] J. Allaby, et al, Zeit. Phys. C36, 611 (1985) . 

[7] E.A. Paschos and L. Wolfenstein, Phys. Rev. D7, 91 (1973) . 

[8] A.J. Buras and K.J.F. Gaemers, Nucl. Phys. B132, 249 (1978) . 

[9] W.G. Seligman, et al, Phys. Rev. Lett. 79, 1213 (1997) . 
[10] M. Arneodo, et al. , Nucl. Phys. B487, 3 (1997) . 
[11] E.A. Hawker, et al. , Phys. Rev. Lett. 80, 3715 (1998) . 
[12] S.A. Rabinowitz, et al. , Phys. Rev. Lett. 70, 134 (1993) . 
[13] CTEQ Collaboration, Phys. Rev. D55, 1280 (1997) . 
[14] L.W. Whitlow, SLAC-Report-357, 109 (1990) . 

[15] D.Yu. Bardin, V.A. Dokuchaeva, JINR-E2-86-260 (1986); and private 
communication. 

[16] J. Pumplin, Phys. Rev. Lett. 64, 2751 (1990) . S < 2 GeV 2 is allowed by 
data summarized in M. Virchaux and A. Milsztajn, Phys. Lett. B274, 
221 (1992) . 



